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Non-Linear  Control  Allocation  Using  Piecewise  Linear 
Functions:  A  Linear  Programming  Approach 

Michael  A.  Bolender  * 

David  B.  Doman  t 

Air  Force  Research  Laboratory,  Wright-Patterson  AFB,  OH  45433  * 


Abstract 

The  performance  of  two  different  approaches  to  solving  the  non-linear  control  allo¬ 
cation  problem  is  presented.  The  non-linear  control  allocation  problem  is  formulated 
using  piecewise  linear  functions  to  approximate  the  control  moments  produced  by  a  set 
of  control  effectors.  When  the  control  allocation  problem  is  formulated  as  a  piecewise 
linear  program,  an  additional  set  of  constraints  enter  into  the  problem  formulation. 
One  approach  is  to  introduce  a  set  of  binary  variables  to  enforce  these  constraints. 
The  result  is  a  mixed-integer  linear  programming  problem  that  can  be  solved  using 
any  branch-and-bound  software.  A  second  approach  is  to  solve  the  piecewise  linear 
programming  problem  using  a  modified  simplex  method.  The  simplex  algorithm  is 
modified  to  enforce  a  subset  of  the  decision  variables  to  enter  into  the  basis  only  if 
certain  conditions  are  met.  We  will  show  that  solving  the  optimization  problem  using 
the  simplex  based  approach  is  significantly  faster  than  solving  the  same  problem  using 
a  mixed-integer  formulation.  We  will  then  compare  the  closed-loop  performance  of  a 
re-entry  vehicle  using  both  approaches. 

^Visiting  Scientist.  Senior  Member  AIAA. 

^Senior  Aerospace  Engineer.  Senior  Member  AIAA. 

*U.S.  Government  work  not  protected  by  U.S.  copyright 


I 


ASC  -  0  4-  01  9  4 


Introduction 


Historically,  the  control  allocation  problem  has  been  approached  by  assuming  that  aerody¬ 
namic  control  effectors  produce  control  moments  that  are  linear  functions  (at  least  locally) 
of  the  control  surface  deflection.  While  this  approach  may  be  valid  for  nearly  all  aircraft 
under  nominal  flight  conditions  when  there  are  no  failed  control  surfaces,  the  occurrence  of 
one  or  more  failed  control  effectors  may  result  in  a  situation  where  the  linear  assumption  is 
no  longer  valid.  When  a  control  effector  fails,  the  resulting  rolling,  pitching,  and  yawing  mo¬ 
ments  produced  by  the  failed  effector  must  be  accommodated  by  reconfiguring  the  remaining 
healthy  control  effectors.  For  some  failures,  one  or  more  of  the  remaining  effectors  may  be 
driven  to  a  position  limit  in  order  to  compensate  for  the  failed  effector.  The  behavior  of  the 
control  moment  curve  in  regions  near  position  limits  tends  be  highly  non-linear.  Therefore, 
the  linear  model  may  introduce  significant  modelling  errors  that  result  in  incorrect  control 
surface  deflections  and  this  may  ultimately  lead  to  loss  of  the  vehicle. 

There  have  been  recent  efforts  to  improve  upon  the  current  control  allocation  approaches 
that  assume  a  linear  relationship  between  the  control  moments  and  the  control  effector.  Do- 
man  and  Oppenheimer1  have  implemented  a  linear  control  allocator  that  added  an  intercept 
term  to  an  estimate  of  the  local  slope  of  the  control  moment  curve  in  order  to  more  ac¬ 
curately  account  for  the  non-linear  behavior  of  aerodynamic  control  effectors.  They  have 
shown  an  improvement  in  tracking  performance  when  there  are  failures  present  without  any 
increase  in  computation  time.  However,  the  intercept  correction  method  assumes  that  the 
control  moment  curves  are  monotonic  functions  of  position,  otherwise  it  is  possible  that 
the  control  allocator  may  settle  in  a  region  where  the  slope  changes  sign,  resulting  in  gross 
modelling  errors  due  to  an  incorrect  intercept  term. 

Doman  and  Sparks2  provided  a  method  for  the  determination  of  the  non-linear  attainable 
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moment  set  (AMS)  for  the  two  moment  case.  Their  effort  considered  control  effectors  that 
produced  a  control  moment  that  was  a  quadratic  function  of  effector  position  about  one  axis 
and  a  linear  function  about  a  second  axis.  More  recently,  Bolender  and  Doman3  extended  the 
work  in  Reference2  to  three  dimensions  for  the  computation  of  the  non-linear  AMS  volume 
when  the  control  moments  about  the  third  control  axis  were  linear  functions  of  effector 
position.  At  the  present  time,  this  particular  methods  does  not  lend  itself  to  be  applied  in 
a  real-time  control  system  due  to  the  excessive  computational  burden  required. 

Most  recently,  Bolender  and  Doman4  have  shown  that  a  piecewise  linear  representation 
of  the  control  moment  curve  accounts  for  the  non-linearities  inherent  in  aerodynamic  data. 
The  only  requirement  is  that  the  control  moments  generated  by  each  effector  be  separable, 
meaning  that  no  aerodynamic  interference  occurs  among  the  effectors.  The  resulting  control 
allocation  problem  was  posed  as  a  mixed-integer  linear  program  and  solved  using  a  public 
doman  branch-and-bound  optimization  code.  The  downside  of  this  approach  was  the  long 
time  needed  by  the  solver  to  find  a  solution,  making  implementation  in  a  digital  flight  control 
system  impracticable. 

The  intent  of  this  paper  is  to  show  that  the  piecewise  linear  control  allocation  approach 
can  be  solved  fast  enough  to  be  implemented  in  a  digital  flight-control  system.  The  approach 
taken  differs  from  the  authors’  initial  paper  on  this  subject  in  two  ways.  The  first  difference  is 
to  move  away  from  the  mixed-integer  linear  programming  form  of  the  optimization  problem 
and  back  to  a  linear  programming  problem  that  imposes  a  set  of  restricted  basis  entry 
rules.  The  second  difference  is  to  solve  a  “mixed  optimization”  problem  as  opposed  to  the 
solution  of  the  multi-branch  control  allocation  problem.  This  allows  us  to  achieve  the  same 
objectives  while  only  having  to  solve  one  optimization  problem  instead  of  two.  We  compare 
the  performance  of  the  simplex  method  with  restricted  basis  entry  rules  to  the  mixed-integer 
formulation  and  show  that  the  two  approaches  give  the  same  performance. 
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Piecewise  Linear  Mixed  Optimization  Control 

Allocation 


The  control  allocation  problem  solved  in  Reference4  was  a  Buffington’s5  multi-branch  ap¬ 
proach.  However,  instead  of  a  linear  relationship  between  the  control  effector  displacements 
and  the  control  moments,  a  non-linear  relationship  was  used.  Let  the  non-linear  vector¬ 
valued  function  G(8)  denote  the  relationship  between  the  control  effectors  and  their  mo¬ 
ments.  The  function  G(8 )  maps  Kn  into  Rm  where  n>  m.  Let  d,ics  denote  the  controlled 
variables.  In  this  case  the  controlled  variables  are  the  rolling,  pitching  and  yawing  moments. 
Buffington’s  multi-branch  approach  requires  that  two  optimization  problems  be  solved.  The 
first  optimization  problem  is  called  the  control  deficiency  branch.  The  objective  of  the  con¬ 
trol  deficiency  branch  is  to  minimize  (C?(5)  —  dde$)\\i  subject  to  position  constraints  on 
<5.  The  value  of  the  performance  index  for  this  optimization  problem  indicates  whether  or 
not  ddes  is  feasible.  If  feasibility  of  the  control  allocation  problem  has  been  ascertained,  a 
second  optimization  problem  is  then  solved  that  minimizes  some  secondary  objective.  The 
objective  function  of  this  second  optimization  is  typically  taken  to  be  ||  w^(8  —  <5p)||i  subject 
to  G(8)  =  ddes  and  position  constraints  on  8.  This  is  commonly  referred  to  as  the  control 
sufficiency  branch. 

The  mixed  optimization  problem  that  was  formulated  by  Bodson6  combines  the  two 
branches  of  the  multi-branch  control  allocation  problem  into  a  single  optimization  problem. 
A  new  parameter  e  is  introduced  for  the  purpose  of  prioritizing  either  control  deficiency  or 
control  sufficiency.  The  mixed  optimization  problem  is  stated  as: 

min  J  =||t»T(G(«)  -  d*s)||i  +  £||«£(i  -  <5^)11,  (1) 

subject  to: 

d min  —  ^  —  8ma x  (2) 
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where  wa  is  an  m  x  1  vector  of  weights  to  prioritize  a  given  control  axis,  G(6)  is  a  non-linear, 
vector- valued  function  that  maps  Rn  — >  Rm,  S  is  an  n  x  1  vector  of  control  effectors,  wu  is 
annxl  vector  of  weights  on  the  control  effectors,  and  Sp  is  an  n  x  1  vector  of  “preferred” 
control  effector  displacements.  Again  it  is  assumed  that  n>m.  For  the  time  being  we  will 
make  the  assumption  that  G(d)  =  BS.  The  mixed  optimization  problem  can  then  be  posed 
as  a  linear  programming  problem.  The  corresponding  linear  program  can  then  be  solved  by 
any  readily-available  linear  programming  software. 


Bodson6  gives  one  possible  transformation  to  a  linear  programming  problem  for  the 
optimization  problem  defined  in  Equations  1  and  2;  however,  we  selected  a  transformation 
approach  that  follows  that  found  in  Bertsimas.7  The  transformation  relies  on  the  observation 
that  |z|  is  the  smallest  number  xs  that  satisfies  both  x  <xs  and  —x  <  xs.  As  a  result,  we 
are  able  to  pose  the  mixed  optimization  problem  as  follows: 

min  J  =  w^6s  +  ew^us  (3) 


subject  to: 


<5,  >0 

(4) 

us  >  0 

(5) 

BS  S3  ^ 

(6) 

—BS  +  Ss  >  “ ddes 

(7) 

S  H-  us  ^  Sp 

(8) 

<5  ^  Sp 

(9) 

S-min  ^  S  <  Smax 

(10) 

The  vectors  6S  and  us  are  vectors  of  “slack  variables” .  The  reason  for  selecting  this  particular 
transformation  as  opposed  to  the  one  in  Bodson6  is  that  this  formulation  allows  us  to  easily 
implement  the  piecewise  linear  function  approximation. 
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To  convert  the  linear  programming  problem  into  a  piecewise  linear  programming  problem, 
we  simply  replace  B5  above  with  a  piecewise  linear  representation  to  each  control  moment 
curve,  (i.e.,  Li(5i),  and  Ni(5i)).  We  choose  a  set  of  breakpoints  such  that  for  each 


*-X3W  (id 

k= 1 

f>S‘>  =  l  (12) 

k= 1 

A,0)Af }  =  0,  ifk>j  +  l,j  =  l,...,Ki-l  (13) 

where  A^  is  a  non-negative  interpolating  coefficient  corresponding  to  the  ith  control  effector 
at  breakpoint  k,  and  Ki  denotes  the  number  of  breakpoints  for  the  ith  control  effector. 
Equation  13  is  necessary  in  order  ensure  that  Si  is  approximated  by  no  more  than  two 
adjacent  values  of  x[k\  If  Si  falls  at  a  breakpoint,  there  will  only  be  one  value  of  A^  that 
is  non-zero  and  Equation  13  is  still  valid. 


The  piecewise  linear  approximations  for  the  control  moments  as  a  function  of  are 
written  as 


(14) 

k=  1 

M,  Mi» 

(15) 

fc=  1 

k=l 

(16) 

where  iff* ,  M\k\  and  Nj®  are  the  values  of  the  rolling,  pitching,  and  yawing  moment  curves 
evaluated  at  the  kth  breakpoint  for  the  ith  control  effector.  We  are  now  able  to  re-write  the 


B  matrix  as 


\l? 

L(!2)  . 

..  . 

T{Km) 

•  •  Tvm 

B  = 

M !(1) 

mJ2)  . 

..  Mlk)  . 

..  Mm™* 

N 

n[2)  . 

..  Nlk)  . 

l \r(Km) 

►  •  d\f}i 

(17) 
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Furthermore,  we  define  a  vector  A  as 


K’l 

A?> 


(18) 


such  that  BS  in  Equations  6  and  7  is  replaced  by  BA.  The  vector  A  is  of  length  ^"_1  A* 

and  B  is  a  matrix  of  size  x  Y!>u=i  K*  where  nCr  is  the  number  of  controlled  variables. 

In  the  piecewise  linear  optimization  problem,  the  constraints  Smin  <  6  <  Smax  are  replaced 

by  \\k)  >  0.  The  upper  and  lower  bounds  on  S  are  accounted  for  in  the  selection  of  the 

* 

breakpoints  for  each  Si.  Once  we  obtain  an  optimal  solution  to  the  problem,  we  compute 
each  Si  using  Equation  11.  It  is  also  necessary  to  include  in  the  problem  the  n  constraints 
that  correspond  to  Equation  12. 


The  resulting  optimization  problem  is 

min  J  =  w^Ss  +  ew^us 


(19) 
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subject  to: 


Ss>  0 

(20) 

us  >  0 

(21) 

BA  +  Ss  >  ddes 

(22) 

—BA  +  Ss  >  —ddes 

(23) 

Ki 

(24) 

t=l 

Ki 

y  y$ik) + us,i  >  -6p>i 

(25) 

i=l 

^ ' 

IV 

o 

(26) 

Ki 

(27) 

1=1 

At0)Af }  =  0,  if  k  >  j  +  1,  j  —  1, Ki  —  1 

(28) 

The  constraint  given  by  Equation  28  forces  us  to  choose  one  of  two  approaches  in  order 
to  obtain  a  solution  to  the  above  optimization  problem.  The  first  approach  is  to  define  a 
set  of  binary  variables  along  with  a  set  of  additional  constraints  that  enforce  Equation  28. 
The  resulting  optimization  problem  is  then  a  mixed-integer  linear  program.  The  second 
approach  is  to  solve  the  optimization  problem  as  a  linear  programming  problem  using  a 
modified  simplex  algorithm.  In  order  to  accommodate  Equation  28,  “restricted  basis  entry 
rules”  are  used  to  determine  which  are  allowed  to  enter  the  basis. 


Mixed  Integer  Linear  Program  Formulation 


Begin  by  considering  the  piecewise  linear  approximation  shown  in  Figure  1.  Note  that  if 
there  are  K  breakpoints,  then  there  are  K  —  1  linear  segments.  We  assign  a  variable 
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that  corresponds  to  the  kth  linear  segment  of  the  piecewise  linear  approximation  such  that 

(29) 

for  k  =  1, . . .  ,K  —  1.  Next,  we  make  the  observation  that  if  A^  ^  0  and  A^  ^  0,  then 


if  A<*>  ^  0  and  A(*+1>  ^  0, 
otherwise 


A(1)  <  yW  (30) 

A(2)  <  yW  (31) 

where  yW  =  1.  However,  if  we  are  on  the  segment  where  A®  ^  0  and  A^  ^  0,  such  that. 
yW  =  1,  then 

A(2)  <  ?/2>  (32) 

A(3)  <  y(2).  (33) 

If  we  proceed  in  this  manner,  we  observe  that  the  following  restrictions  can  be  placed  on 
the 


A  {1)<y(1\  (34) 

X(k)  <  +  y{k),  k  =  2, . . . ,  K  -  1  (35) 

AW  <  y(K-l)  (36) 

The  rationale  behind  Equation  35  is  as  follows:  the  A^  that  correspond  to  points  that  are 
interior  to  the  interval  (i.e.,  they  are  not  the  endpoints  of  the  interval  on  which  x  is  defined) 
can  be  associated  with  one  of  two  fine  segments.  A  particular  A^  is  the  endpoint  for  the  line 
segment  immediately  preceding  it  in  addition  to  the  line  segment  that  comes  immediately 
after  it.  Only  one  of  these  two  line  segments  may  be  “active”  at  any  time;  therefore,  the 
right-hand  side  of  Equation  35  is  never  greater  than  one.  In  addition  to  Equations  34-36,  we 
have  an  additional  constraint  to  ensure  that  only  one  of  the  K  —  1  line  segments  is  active, 
hence  only  one  of  the  y ^  can  be  equal  to  one: 

X>(fe)  =  1  (37) 

k= 1 
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Figure  1:  Piecewise  Linear  Function  Approximation 
The  above  can  be  generalized  for  the  multi- variable  case  as  follows. 


\-1}  <  y^\  i  =  1, . . .  ,m 

(38) 

\-fe) 1  <  yf ~1}  +  yf‘ \  *  =  l,....,n»,  k  =  2, . . .  ,Ki  —  1 

(39) 

Af*  <  *  =  1, . .. ,m 

(40) 

Ki- 1 

E^_1)  =  1.  <  = 

(41) 

A;=l 

V?: ’  €  {0, 1} 

(42) 

Equations  38-42  replace  Equation  28  above.  As  a  result,  we  have  a  mixed-integer  linear 
programming  problem. 


Linear  Program  with  Restricted  Basis  Entry  Rules 

1 

An  alternate  method  of  handling  Equation  28  is  to  solve  the  piecewise  linear  program  given 
in  Equations  19-28  minus  Equation  28.  Without  Equation  28,  we  are  left  with  a  linear 
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programming  problem  that  can  be  easily  solved  using  the  simplex  method.  However,  we  need 
an  alternate  approach  of  enforcing  the  adjacency  constraint.  Instead  of  defining  Equation  28 
explicitly  and  including  it  in  the  problem  formulation,  it  is  to  be  enforced  through  a  set  of 
rules  in  the  simplex  algorithm  that  only  admit  a  if  Equation  28  is  satisfied.  The  rules 
are  commonly  referred  to  as  “restricted  basis  entry  rules”.  More  details  on  the  simplex 
algorithm  with  restricted  basis  entry  rules  can  be  found  in  Miller8  and  Hadley.9 

A  survey  of  commercial  linear  programming  solvers  indicated  that  simplex  based  solvers 
with  restricted  basis  entry  rules  were  unavailable.  The  GNU  Linear  Programming  Kit 
(GLPK),  which  is  an  open-source  general  linear  program  and  mixed-integer  linear  program 
solver  written  in  C,  was  modified  to  accommodate  restricted  basis  entry  rules. 


Results 


The  results  presented  in  this  section  are  given  for  an  unpowered  re-entry  vehicle.  This 
particular  vehicle  utilizes  the  following  six  control  effectors:  left  and  right  ruddervators,  left 
and  right  flaperons,  a  speed  brake,  and  a  body  flap.  We  show  that  the  mixed-optimization 
problem  implemented  with  piecewise  linear  approximations  of  the  control  moments  and 
solved  using  a  simplex  method  with  restricted  basis  entry  rules  gives  a  solution  that  is 
comparable  to  that  given  in  Bolender  and  Doman.4  We  also  demonstrate  that  the  mixed- 
optimization  problem  can  be  solved  fast  enough  that  it  is  a  candidate  for  future  use  on  a 
digital  flight  control  computer  equipped  with  a  processor  that  is  comparable  to  one  found  on 
current  desktop  computers.  For  this  study,  the  GLPK  solver  was  implemented  as  a  C  mex 
file  and  compiled  to  minimize  execution  speed.  The  computer  on  which  this  analysis  was  run 
was  equipped  with  an  Athlon  1800XP+  processor,  1.5  GB  of  RAM,  and  the  Windows  2000 
operating  system.  We  compare  the  execution  time  of  the  simplex  method  with  restricted 
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basis  entry  rules  to  the  time  required  to  solve  the  same  optimization  problem  using  the 
mixed-integer  linear  programming  formulation  of  the  mixed-optimization  problem.  The 
mixed-integer  linear  programs  are  solved  using  the  branch-and-bound  solver  supplied  by  the 
GLPK  package. 


Simulation  Results 

The  results  that  follow  give  the  closed-loop  vehicle  performance  when  there  are  two  failures 
injected  into  the  flight  control  system  at  different  times  during  the  approach  and  landing 
phases.  It  is  assumed  that  there  is  some  type  of  fault  detection  capability  on-board  the 
aircraft  to  identify  the  failures.  The  failure  information  is  immediately  passed  to  the  control 
allocation  algorithm  in  order  to  facilitate  re-configuration  of  the  vehicle’s  effectors.  The 
aircraft’s  trajectory  begins  at  an  altitude  of  about  15,000  ft  above  the  runway  and  4  miles 
downrange  from  the  runway  threshold.  The  first  failure  occurs  30  s  into  the  simulation,  and 
involves  the  body  flap  being  locked  at  —5  deg.  This  failure  contributes  a  constant  pitching 
moment  to  the  aircraft.  A  second  failure,  where  the  right  rudder  becomes  locked  at  1  deg, 
occurs  at  40  s.  This  particular  failure  adds  not  only  a  pitching  moment  to  the  aircraft,  but 
also  rolling  and  yawing  moments.  This  particular  failure  combination  was  chosen  because 
it  requires  the  flaperons  to  operate  in  a  highly  non-linear  region  of  the  control  moment 
curve.  After  the  failures  are  introduced,  the  aircraft  tries  to  follow  the  nominal  approach 
trajectory  to  the  runway  threshold.  The  aircraft  extends  the  landing  gear  at  about  68  s 
and  flares  immediately  before  touchdown.  The  simulation  ends  at  touchdown  when  the 
weight-on-wheels  switch  is  triggered. 

Presented  below  in  Figures  2  and  3  are  the  time  histories  for  the  control  moment  error 
and  the  control  effector  command  time  histories.  In  Figure  2  the  moment  error  is  defined  by 
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log10  || BA  —  ddes\[2  where  B  and  A  are  defined  above.  Recall  that  B  is  the  piecewise  linear 
analog  to  the  B  matrix  that  one  encounters  when  using  a  linear  approximation;  therefore, 
the  product  BA  is  the  linear  interpolation  of  the  control  moment  data.  We  see  that  the 
modelling  errors  are  negligible  until  the  rudder  failure  is  introduced  at  the  40s  mark.  Between 
40s  and  60s  there  is  a  control  deficiency;  therefore,  the  desired  moment,  djes  is  not  feasible 
under  the  failure  conditions.  The  performance  of  the  mixed-optimization  is  identical  to  that 
given  in  Reference  [4]. 

The  control  surface  time  histories  given  in  Figure  3  are  the  control  deflections  returned 
by  the  control  allocator.  The  control  surface  deflections  compare  favorably.  There  are  very 
minor  differences  in  the  right  flap  deflection  after  the  rudder  failure  is  introduced,  but  this 
discrepancy  does  not  appear  to  be  an  issue. 


Performance  Results 

The  result  given  in  Table  1  is  for  a  fixed  flight  condition,  and  considers  10000  different  control 
input  vectors  with  each  component  selected  randomly  from  a  uniform  distribution.  For  each 
input  vector,  the  control  allocation  problem  was  solved  using  both  a  mixed-integer  linear 
program  and  simplex  with  restricted  basis  entry  rules.  The  preference  vector  was  taken  to 
be  Sp  —  0.  The  execution  time  of  both  approaches  was  measured  using  the  standard  C 
library  function  clock ().  It  is  apparent  that  the  simplex  method,  on  average,  is  an  order 
of  magnitude  faster  than  the  MILP  approach.  If  it  is  assumed  that  the  typical  sample  time 
of  a  digital  flight  control  system  is  0.02  sec,  then  the  solution  time  given  in  Table  1  for  the 
simplex  method  with  restricted  basis  entry  is  more  than  adequate  for  practical  application 
given  a  processor  that  is  in  the  same  class  as  that  tested  above.  The  mean  control  deflection 
for  each  method  is  given  in  Table  2.  Note  that  for  this  particular  set  of  random  ddea,  the 
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Figure  2:  Difference  Between  Commanded  and  Applied  Moments 
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Left  Rudder 


Surface  Deflections 


Mean 

Std.  Dev. 

MILP 

0.0885 

0.0423 

RBE 

0.0083 

0.0042 

Table  1:  Solver  Execution  Time  Statistics:  Athlon  XP1 800+/ Windows  2000 


Mean 

Std.  Dev. 

MILP 

18.2365 

8.6293 

RBE 

18.2657 

8.5736 

Table  2:  Mean  and  Standard  Deviation  of  ||d||2. 


mean  deflections  are  nearly  the  same  for  both  approaches.  On  the  other  hand,  Table  3  shows 
a  significant  difference  in  the  average  moment  error.  There  were  a  small  number  of  d^s  for 
which  the  Simplex  Method  with  Restricted  Basis  Entry  Rules  failed  to  find  a  set  of  control 
deflections  that  would  produce  a  feasible  moment  while  the  MILP  formulation  succeeded. 
The  average  error  for  these  vectors  is  0.4938  and  is  the  cause  for  the  poor  correlation  between 
the  MILP  solution  and  the  Simplex  with  Restricted  Basis  Entry  Rules.  For  each  of  these 
vectors,  the  control  deficiency  occurs  in  the  yawing  moment.  Although  it  occurs  in  less  than 
0.5%  of  the  cases  that  were  tested,  the  cause  of  this  behavior  is  currently  under  investigation. 


Conclusions 


The  approach  to  non-linear  control  allocation  that  was  presented  assumed  that  control 
moments  generated  by  the  deflection  of  aerodynamic  surfaces  were  separable  functions. 
This  assumption  allows  us  to  approximate  any  separable  non-linear  function  by  a  piecewise 


Mean 

Std.  Dev. 

MILP 

4.4547  x  10"^ 

1.0709  x  10~lbr 

RBE 

1.4139  x  10"5 

2.6927  x  10~4 

Table  3:  Mean  and  Standard  Deviation  of  ||£?A  —  ddes lb- 
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linear  function  in  the  control  allocation  problem.  As  a  result,  we  were  able  to  cast  a  non¬ 
linear  optimization  problem  as  a  linear  programming  problem  where  a  subset  of  the  decision 
variables  are  subject  to  a  set  of  restricted  basis  entry  rules.  Although  this  still  gives  an 
approximate  solution  to  the  control  allocation  problem,  it  is  much  more  accurate  than  the 
traditional  methods  that  assume  linear  relationships  between  the  control  moments  and  the 
control  effector  positions  .  It  was  subsequently  shown  that  a  simplex  algorithm  that  enforces 
the  restricted  basis  entry  rules  is  probably  fast  enough,  given  the  appropriate  processor  clock 
speed,  for  use  in  a  real  time  flight  control  system. 
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